't Hooft Loops, Electric Flux Sectors and Confinement in SU{2) Yang-Mills Theory 
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We use 't Hooft loops of maximal size on finite lattices to calculate the free energy in the sectors 
of SU{2) Yang-Mills theory with fixed electric flux as a function of temperature and (spatial) vol- 
ume. Our results provide evidence for the mass gap. The conflnement of electric fluxes in the low 
temperature phase and their condensation in the high temperature phase are demonstrated. In a 
surprisingly large scaling window around criticality, the transition is quantitatively well described 
by universal exponents and amplitude ratios relating the properties of the two phases. 



PACS numbers: 12.38.Gc, 12.38.Aw, 11. 15. Ha 

Center symmetry is widely believed to play a key role 
for confinement. In SU{N) Yang-Mills theory at finite 
temperature, the Polyakov loop is commonly used in lat- 
tice studies to illustrate this role. Its correlations are 
short range at low temperature and acquire a non-zero 
disconnected part above a critical temperature Tc, sig- 
naling deconfinement. The role of the low (high) tem- 
perature phase as the center symmetric (broken) one is 
thereby the opposite from that in the corresponding Z^r- 
spin model. This suggests to consider dual variables, 
whose behavior as a function of temperature is reversed. 
That point of view was emphasized by 't Hooft with the 
introduction of the dual to the Wilson loop ||l| . Temporal 
't Hooft loops in SU{2) show screening for the interac- 
tion of a static pair of center monopoles in both phases, 
for T <Tc and for T > T^ [|,|, just as spatial Wilson 
loops exhibit an area law in either case. The expectation 
values of (sufficiently large) spatial 't Hooft loops W{C), 
on the other hand, change from a screening behavior be- 
low Tc to a confined one 0] with a dual string tension 
a and an area law in the "electrically" deconfined phase 
above T, @, 

{W{C)) ~ exp{-CT(r) LR} , at T>Tc, (1) 

for a rectangular curve C spanning a spatial surface of 
size L X R on a. 1/T X L^ lattice. Ref. [pi also confirmed 
numerically a perimeter law for 't Hooft loops at T = 0. 

The qualitative behavior of the spatial 't Hooft loops 
is thus the same as that of the Wilson loops in the 3- 
dimensional Z2 gauge theory Q. Furthermore, as the 
phase transition is approached (from above), the tem- 
perature dependence of the SU{2) dual string tension 
obeys the same scaling law as that of the interface ten- 
sion in the 3-dimensional Ising model (below Tc). This 
similarity is expected, since the 3d Z2 gauge theory, its 
dual the 3-dimensional Ising model, and SU{2) at finite 
temperature all belong to the same universality class. 

To introduce a spatial 't Hooft loop of maximal size 
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L X L, living in, say, the {x,y) plane of the dual lat- 
tice, one multiplies one plaquette in every {z,t) plane of 
the original lattice by a non-trivial element of the cen- 
ter of SU{N) in such a way that the modified plaquettes 
form a coclosed set. This creates a Zjv-interface which is 
equivalent to enforcing boundary conditions with twist in 
the (z, t) directions. Combining two and three maximal 
LxL\ Hooft loops in orthogonal spatial planes yields the 
partition functions of SU{N) Yang-Mills theory on finite 
lattices for all possible combinations of temporal twists. 
From these, one obtains the free energies in the presence 
of fixed units of electric fiux e € Z^ via a Zjv Fourier 
transform as shown by 't Hooft |^ . This leads to the re- 
spective electric-flux superselection sectors of the theory 
in the thermodynamic limit. In the present paper, we 
verify for SU{2) that the partition functions of electric 
fluxes e ^ vanish exponentially with the spatial size L 
of the system for T < Tc in the magnetic Higgs phase 
with electric conflnement, whereas for T > Tc they be- 
come equal to that of the neutral e = sector. Describing 
the transition by critical exponents of the 3d-Ising class, 
we discuss how universal amplitude ratios quantitatively 
relate the two phases, in particular, the string tension 
below and the dual (vortex) string tension above Tc- 

As expected, the free energy Fg of a static fundamental 
charge jumps from -l-oo to at Tc- One might expect to 
see this also in measuring the Polyakov loop P directly 
on the lattice. If (P) = e~T^'', an inflnite free energy 
amounts to a center symmetric distribution, while a non- 
zero expectation value (P) is obtained for flnite Fq. How- 
ever, the presence of a single charge is incompatible with 
periodic boundary conditions to measure (P). And, like 
any Wilson loop, (P) is subject to UV-divergent perime- 
ter terms, such that (P) = at all T as the lattice spacing 
a—fO. Here, we measure the gauge- invariant, UV-regular 
free energy of a static fundamental charge in SU{2), and 
show that it has the expected behaviour, dual to that of a 
certain type of center vortex. Both provide a well-defined 
order parameter for the transition |^]. 



Twisted Boundary Conditions and Electric Fluxes 

For the finite-volume partition functions of the pure 
SU{N) gauge theory, 't Hooft's twisted boundary con- 
ditions fix the total number of ZTv-vortices modulo N 
that pierce planes of a given orientation. Thus, on the 
4-dimensional torus there are N^ different Z^r-fiux sec- 
tors corresponding to the 6 possible orientations for the 
planes of the twists. These label the inequivalent choices 
for imposing boundary conditions on the gauge potentials 
A, which are invariant under the center Zat of SU{N). 
One first chooses A{x) to be periodic with the lengths 
Lf^ of the system in each direction /x up to gauge trans- 
formations ri^(xx) G SU{N) which can depend on the 
components x±_ transverse to that direction, 

A{x+L^) = n^{x^) {A{x) ~-d) n-\x^) . (2) 

Then, compatibility of two successive translations in a 
(/i, i^)-plane entails that (no summation of indices) 



i^f,{x±+L^)Q^{xi_) = Zf^^D,^{x±+Lfj^)Qf,{x±_) 



with Zf,^ = e^^^n^JN ^ ^^^ ^^^ 



-N 



(3) 



The total number mod. N of center vortices in a (/i, ly)- 
plane is specified in each sector by the corresponding 
component of the twist tensor n^i, . The spatial ones are 
given by the conserved, Z^v-valued and gauge-invariant 
magnetic flux m through the box, Uij = tijumu- The 
time components n^i = ki define temporal twist k S Z^. 
With the inequivalent choices of boundary conditions, 
the finite- volume theory decomposes into sectors of frac- 
tional Chcrn-Simons number (t^ + A; • rn/N) |^ and states 
labelled by |fc, m, i/), where i^ e Z is the usual instanton 
winding number. However, these sectors are not invari- 
ant under homotopically non-trivial gauge transforma- 
tions 17 [A:, I/] which can change k and v, 



Q[k',iy']\k,rn,iy) == \k + k' ,rn,iy+iy') 



(4) 



A Fourier transform of the twist sectors Z(k, rn, v) which 
generalizes the construction of 6'-vacua as Bloch waves 
from J/- vacua in two ways, by replacing v —^ (y + k-rfi/N) 
for fractional winding numbers and with an additional 
Z^v-Fourier transform w.r.t. the temporal twist fc, 



g-if(e>i.9) = j^Yl e-'^^^^'"^ Z{k, 



(5) 



yields the free energy -F(e, m, 6) in an ensemble of states 
invariant, up to a geometric phase uj(k, v) = 27re • k/N + 
0{v+k ■ m/N), under the non-trivial Q[k, v] also. 



il[fc, J/] |e, m, 0) — exp{iti;(fc, ;/)} |e, TO, 



(6) 



These states are then classified, in addition to their mag- 
netic flux TO and vacuum angle 0, by their Zjy-valued 
gauge-invariant electric flux in the e-direction [^. 
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FIG. 1. Cubes with one, two, and three L x L planes dual 
to the stacks of flipped plaquettes, sketched for temporal twist 
k = (0, 0, 1), k = (0, 1, 1), and k = (1, 1, 1) from left to right. 

To create n^j,-twist in SU{2), on the lattice, one intro- 
duces couplings with reversed signs into the usual Wilson 
action by replacing /3 — + — /3 for one plaquette in every 
(/i, i^)-plane. 



5(/3,fc,TO) = -^/3(P)iTr(t/p) 



(7) 



Herein, the sum extends over all plaquettes P with Up 
denoting the path-ordered product of the links around P 
and we introduced a plaquette-dependent coupling, 



/3(P) = 



-(3, PeVin^,) 
/3, P^PK.) 



(8) 



where 'P{n^^) denotes the coclosed stacks of plaquettes 
dual to the planes of the maximal 't Hooft loops. These 
stacks of flipped plaquettes force the Z2-interfaces corre- 
sponding to twist in the (/i, j/)-directions. Equivalently, 
they create a 't Hooft loop of maximal size in the orthogo- 
nal plane. For the various combinations of temporal twist 
the coclosed sets 'P(fc), when put between two time-slices, 
can be chosen dual to the spatial planes shown in Fig. ^. 
The partition functions of the twist sectors, relative to 
the untwisted Zp (such that Z/3 ((),()) = 1), are then: 



Zfs{k,rn) — Z„ / \dU\ exp — S'(/3, fc, to) 



(9) 



Similarly, flipping the couplings of plaquettes dual to 
some surface subtended by a closed curvs^^J^ yields the 
expectation value of a finite 't Hooft loop WiC). An en- 
tirely analogous procedure can be used for Wilson loops 
in the 3-dimensional Z2 gauge theory. Through dual- 
ity, their expectation values can be expressed as ratios of 
Ising-model partition functions with and without antifer- 
romagnetic bonds at those links of the Ising-model that 
are dual (in 3 dimensions) to some surface spanned by 
the Z2-Wilson loop g . In both cases the surface is arbi- 
trary except for its boundary. The 't Hooft loop can thus 
be viewed as a gauge-invariant operator which creates a 
fluctuating center-vortex surface with pinned boundary. 
In this paper, we first calculate Z^i^) = Zp{k,0) for 
TO = (and 6* = 0) on a 1/T x L^ lattice in SU{2) 
(with ki in {0, 1}). The expectation values of maximal- 
size 't Hooft loops, given by the partition functions of 



the twist sectors, are then used to calculate the free 
energies of electric fluxes as per Eq. (|5|). For purely 
temporal twists in particular, the free energies of the 
electric fluxes through the L^ box at temperature T, 
Fe(e; L, T) = F(e, m = 0, 6* = 0) - F(e = 0, m = 0, 6* = 0), 
are given by: 



^e(e) 



^ -l,F,(e-L,T) 






'■"/'' Zk{k) 



Effc.io} ^k{k) 



(10) 



with Ze(d) — Zk{0) = 1. Because of the invariance un- 
der spatial tt/2 rotations, we can write for SU{2): Zk{\) 
if fc = {(1,0,0), (0, 1,0), (0,0, 1)}; Zk{2) if fc = {(1, 1,0), 
(1, 0, 1), (0, 1, 1)}; and Zfc(3) = Zfc(l, 1, 1), for the parti- 
tion functions with one, two and three maximal 't Hooft 
loops in orthogonal spatial planes, respectively. With 
analogous notations for the Zp{e) one thus obtains: 



Ze(l) = ;;^( 1 + Zk{l) - Zk{2) - Zfe(3) 



Z,{2) = ;^( 1 - Zk{l) - Zk{2) + Zfe(3) 



Ze(3) = ;;^(1 - 3^fc(l) + iZk{2) - Zk{i) 
with J\fe ^ 1 + iZkil) + 3Zfc(2) + Zfc(3) 



(11) 
(12) 

(13) 



Eqs.([ll|)-(|l3|) are readily inverted via inverse Z2 Fourier 
transform, which in effect interchanges Ze{i) <-^ Zk{i). 

We now establish the connection with Polyakov loops. 
First, recall that the gauge invariant definition of the 
latter in the presence of temporal twists is given by Q , 



P{x) = Itr (Ve'^ 



rl/T 

j^ Ao{x,t)dt 



^tix) 



(14) 



Successively transforming the path-ordered exponential 
by the various spatial transition functions which accom- 
pany the possibly non-trivial ^t{x) for the transition in 
the time direction, we derive from (H) and (g) that 



P{x + Le) = e^""'^-^/^ P{x) , or 
P{x)P\x + Le) = e-2^"--'=/^ ]| . 



(15) 



This is proportional to the unit operator when acting 
on the states in a sector of definite fc-twist. Therefore, 
the electric-flux partition functions of Eq. ( p^ ) are in fact 
the expectation values of Polyakov loop correlators in the 
ensemble average over all these temporal twists. 



^e(e) 



-^F^{e;L,T) 



{P{x)P\x + Le)) 



L,T 



(16) 



This expectation value is taken in the no-flux ensemble. 



N-l 



with enlarged partition function Z — X^u- =0) Zk{k), 
which is manifestly different, in a finite volume, from 
the periodic ensemble. Also note that the operator in 



Eq. (16) has no perimeter, is UV-regular, and we will see 
that there is no Coulomb term for small volumes either. 
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FIG. 2. The free energy of one temporal twist as a function 
of the finite size scaling variable x (with a; < for T < T^). 

Eq. ( [lO| ) thus yields a dual relation between Polyakov 
loops and temporal twists of the general pattern. 



{P{x)P\x + Le)) 



0, forZfe(fc) 

1, forZfe(fc) 



1 , T < T, 
, T > T, 



reflecting the different realizations of the electric Zjv cen- 
ter symmetry in the respective phases. 

Finite Size Scaling and Universal Amplitude Ratios 

We compute the three partition functions Zk{i) for 
SU{2) near Tc by Monte Carlo, using the method of Ref. 
H, with 20k"50k measurements per simulation. For pi- 
oneering related work at T = 0, see Ref. [nOJ. The Zk{i) 
are the analogues of ratios of 3d-Ising model partition 
functions with different boundary conditions. As for the 
latter ||l^, we assume their L, T dependence to be gov- 
erned by simple finite-size scaling laws. 



Zk{i) 



J± 



(^), 



1,..3. 



(17) 



.(■0 



The /j_ are functions of the finite-size scaling variable 

X = ±Lr,|C oc L/e±(i) , for T ^ T, , (18) 

where t = T/Tc - 1 and ^±{t) = £,1^" are the reduced 
temperature and the correlation lengths, respectively, 
and we use the exponent v = 0.63 from the Ising model. 
In our calculations we keep the number of points Nt in 
the time direction fixed and control the temperature by 
varying the lattice coupling P around the critical value 
for the phase transition f3c. The results presented here 
are obtained with Nt ^ 4 for which Pc = 2.29895(10) [||. 
We employ couplings /3 between 2.19 and 2.5 for various 
lattices ranging from A^j = 6 up to Ni = 20 points in the 
spatial directions. Our method to calculate a partition 
function for one temporal twist amounts to N^ indepen- 
dent Monte Carlo simulations |3| . For a fixed statistical 
accuracy, which introduces another factor A'^;^, the cost 
of the calculation roughly increases as Nj . 



For the temperature T = l/{Nta), where the lattice 
spacing a = a{(3) depends on the coupUng, we adopt the 
leading scaling behavior around criticality of the form, 



T/T, = exp{& (fi - /3e)} 



(19) 



The non-perturbative coefficient herein, we use b — 3.26, 
is determined so as to reproduce the published values of 
Pc{Nt) also for iVt = 6 and 8 ||. 

We can see in Fig. that all our curves intersect at the 
same point T = Tc, which is consistent with the known 
value of /3c, even for our smallest lattice size. In fact, this 
intersection point is known to provide a quite accurate 
determination of the critical coupling f3c already on very 
small lattices in the Ising model [|ll| . 

Away from criticality, corrections to ( jH^ ) will become 
important if one is interested in a more precise definition 
of the physical temperature for SU{2), as in Ref. [[l3[ . 
Especially at our lowest (3 values, where such corrections 
are noticable, scaling violations become important also. 
A more refined definition of the reduced temperature t is 
therefore beyond the scope of our scaling analysis which 
concerns the dominant behavior near Tc- It has no effect 
on the qualitative conclusions emphasized here. 

That our results for all different lattice sizes nicely col- 
lapse on a single curve can be seen for one spatial 't Hooft 
loop in Fig. g, with analogous results for i = 2 and 3. We 
fit the free energies by an Ansatz for each phase, where 
the two leading terms obey the expected thermodynamic 
behavior, while the others represent an ad-hoc modeling 
of small-size corrections (with a^' ~ c[^ < 6^'; d^' < 1), 



Ft\-) = 



cxp 



{6«:r 



,W 



6^'a;2 



,W 



+ 



X <0 



(20) 



x>0 



and 



ln/«(x) 



F. 



^'\x) for a; ^ 0. 



Plotted over tem- 



perature, the data of Fig. || and the unique function 
F^^'(a;) lead to the family of curves shown in Fig. g (top) 
in which the phase transition is exhibited most clearly. 
The amplitudes relevant to the large size behavior of the 
free energy near criticality come out as 6*:^' = 3.87 ± 0.5 
and &^' — i7q^' = 5.36 ±0.1, with some additional system- 
atic uncertainty inherent in the form of our Ansatz (20). 
With analogous data and fits for Zk{2) and Zk{3), from 
Eq. (11), we obtain Ze{l) as shown in Fig. H (bottom). 
Corrections to scaling do not become appreciable up to 
T ^ 2Tc indicating a surprisingly large scaling window. 
Above Tc, the dual string tension is (with CTq^' = 6^'), 



a(r) 



\tr = R/e+it) 



(21) 



where the universal ratio R ~ 0.104 [|lj,|l5[ is known 
from the 3d-Ising model. There, R = £,^ai relates the 
correlation length and the interface tension aj for T < Tc- 
Here, Eq. ( pi| ) determines the screening length for the 



Polyakov loops above Tc, £,+ {t) = yjR/a{T). In addi- 
tion, the universality hypothesis relates the ratio of the 
correlation lengths for the Polyakov loops in SU{2) to 
that of their dual analogue, the correlation lengths of the 
spins in the 3d- Ising model, as measured in Ref. |15||, 



>S(7(2) 



,^SU(2} 



er'/er"' 



1.96 



(22) 



Together with (Ell) this relates the string tension ampli- 
tude below to its dual counter part above Tc, as follows. 
From the linear part of the electric-fiux free energy, 

-<t(T)L/T ^g-L/C_(t)^ T<Tc 

-x/(d"^rc) , - 



{P{x)PHx + Le,)) 

^ -HZe{l)) ^ -x/iC'Tc) , 'X = LTcir, (23) 
for large L (or (—a;) large). Thus, from Eqs. (pl|), ( p2[ ) 
aiT) 1 



T 



?-W 



Tr 



lirv^^^. 



i?4 



e 

^R. 



0.4 



A value of about o-q' ~ 5.36 then implies for the string- 
tension amplitude 1/(^1l'Tc) ~ 3.66. This value was used 
in the fit for the slope of the linear part of — ln(Ze(l)) 
shown in Fig. |j. Fitting the slope to the data yields the 
consistent value l/{£,''^^Tc) = 3.58 ±0.5. This is the linear 
potential between static charges without Coulomb part. 
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FIG. 3. The partition functions of one temporal twist (top) 
and one electric flux (bottom) over T for the various lattices. 
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FIG. 4. Free energy of one unit of electric flux from Eq. U/U 
as a function of x in the confined phase, and above Tc (insert) . 

The insert of Fig. shows the quahty of the finite size 
scaling from the more accurate data for the free energy 
of one electric flux in the high temperature phase. 

Similar results are obtained from Eqs. (|l^), ( [l3[) also 
for 2 and 3 orthogonal electric fluxes which we verify 
to be suppressed more strongly in the confined phase, 
Zei2)/Zeil},Zj3)/Ze{l) ^ for L ^ oo. Then, invert- 
ing Eqs. ([ll|)-([13|) one therefore deduces 



- ln(Zfc(l)) ^ - ln(l - 2Z,(1)) ^ Z,{1) 



i/(d"¥,: 



= 6*^' =3.87 ±0.5 



T<T, 

(24) 



for the string tension, which is again consistent with the 
value ~ 3.66 implied by the universal amplitude ratio. 
Due to similar relations for Zfc(2) and ^^(3), we expect 
6^' = 6^' = 6^^' in our fits (|o|) for 2 and 3 spatial 't Hooft 
loops below Tc, which is verified well within our 0(10%) 
accuracy on these amplitudes. Above Tc on the other 
hand, we expect for the dual string tension amplitudes 



;^(l) . A2) . ~(3) 



C, 







0-, 











1 : V2: \/3 



(25) 



according to the effective area of diagonal loops. In Fig. H 
such a square-root behavior is successfully enforced on 
the fits (^ of the spatial 't Hooft loops above Tc- We 
obtain the same ratios, with less accuracy, for the slopes 
of the electric-flux free energies below T^ as expected for 
diagonal fluxes with string formation. 

To summarize, we have shown that, below Tc, the free 
energy of electric fluxes diverges linearly with the length 
L of the system. Because spatial twists share their qual- 
itative low-temperature behavior with the temporal ones 
considered here, the free energy of the magnetic fluxes 
must vanish. This is the magnetic Higgs phase with elec- 
tric confinement of SU{2) Yang-Mills theory. 

At criticality all free energies rapidly approach their 
finite L —^ oo limits indicative of massless excitations. 
We obtain, e.g., Zk{l) = 0.54(1) for T = Tc, which agrees 
with the corresponding ratio in the 3d-Ising model |15[. 
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FIG. 5. Free energy of 1, 2 and 3 't Hooft loops at T > Tc. 

Above Tc, the free energy of electric charges vanishes 
in the thermodynamic limit. The dual area law prevents 
magnetic charges from propagating in spatial directions. 

The transition is well described by exponents and am- 
plitude ratios of the 3d-Ising class, see also |16|] . 
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